Propagation of Fronts and Information in Dispersive Media 
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We present a general proof based on Kramers-Kronig relations that, in a normal or anomalous 
dispersive linear medium, any (discontinuitynonanalytic disturbance) in an electromagnetic pulse 
can not propagate faster than the phase velocity, c. Consequently the information carried by the 
discontinuity (nonanalytical disturbance) can not be transmitted superluminally. 
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According to the Einstein's theory of special relativ- 
ity, the speed of any moving object can not exceed the 
speed of light in vacuum c. However it is well-known 
that the group velocity of a light pulse Vg can exceed c 
in an anamolously dispersive medium 0]. This interest- 
ing effect is a result of interference of different frequency 
components of the light pulse 2, 3J. The superluminal 
phenomenon disappears when the pulse loses the coher- 
ence Q. 

Sommerfeld and Brillouin observed in 1914 that a 
superluminal group velocity does not violate causality 
[MlEIZi' They observed that the front velocity, (the veloc- 
ity of a sharp nonanalytic discontinuity in a light pulse), 
should be used as the signal velocity at which informa- 
tion is transmitted and this velocity does not exceed c. 
The analysis of Sommerfeld and Brillouin is based on the 
propagation of an electromagnetic wave of the form 
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through an anamolously dispersive medium with strong 
absorption, characterzed by a Lorentzian susceptibilty 
that is proportional to w^^. The sharp begining of the 
light wave (1) corresponds to the signal. The debate con- 
cerning the information velocity in media still remains 
due to progress in experiments. For example, a recent 
experiment by Wang et al. reports superluminal propa- 
gation in a gain medium with susceptibility proportional 
to uj~^ 1^. It is therefore a matter of great interest to 
give a general proof about the causal nature of the prop- 
agation of classical information (carried by the front or 
discontinuities) with a velocity less than or equal to c. 

We note that the peak superluminal propagation of 
a light pulse in a dielectric medium, when the spec- 
trum of the incident pulse is in the anomalous disper- 
sion frequency range of the medium, can be derived 
from Maxwell's electromagnetic theory. The superlu- 
minal propagation is therefore a classical phenomenon 
and the question whether the information can be trans- 
mitted with a velocity faster than c needs to be ad- 



dressed classically. Any practical pulse must have a be- 
ginning (a star ting point of a non-equilibrium process) 
Furthermore, any nonanalytical dis- 
turbance (a discontinuityof the field or its first or higher 
order derivatives) |^. [loLlT^ carries classical information. 
The front of the pulse is one of such discontinuities. 

In this letter, we prove, based only on Kramers-Kronig 
relations and the Maxwell equations, that any nonanalyt- 
ical disturbance (any discountinuity including the front) 
in pulses propagates at the phase velocity in a linear 
medium. It should be emphsized that our general and 
rigorous proof has no requirement on the form of inci- 
dent field and the particularity of the medium. 

We consider the propagation of an electromagnetic 
pulse through a medium occupying the spapce from z = 
to an arbitrary z > whose response to the electric field 
of the light pulse is linear. After passing through the 
medium, the light field at the position z and at time t 
can be written as 
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where 
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has the property of the retarded Green function. Here 
i?o(0, ti) is the input field and = n(ijj) — 1-1- iniuj) 
with n[ijj) and k[ijj) being the real and imaginary parts 
of the complex refractive index, respectively. Since n[ijj) 
and k{uj) satisfy the Kramers-Kronig relations ISj, it can 
be shown that a{uj) is an analytic function in the upper 
half plane and consequently cr(w) can be expanded into 

oo 

the form: cr{Lj) = ^ An/uj^ for — + cx3 , where An are 

n=l 

the expansion coefficients. 

For the pulse propagation through vacuum, (j{uj) = 0. 
Therefore the output field at z is Ey{z,t) — Eo{0,t — 
2:/c),i.e., the light pulse propagates with velocity c. 
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For the pulse propagation through a medium, the func- 
tion G'(^; z) as given by Eq. (3) can be rewritten as 



G(C;z) = ^ I dLoe^'^^ 
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where 
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We then have 

(z, t) = e'^Eoiz, t - z/c) + e'^ E^iz, t), (6) 



where 



Ejrt{z,t) = 



t-z/c 



dtiEa{Q,ti)J{ti -t + z/c;z) (7) 



From Eq. we see that the output field has two parts: 
the first part is an instant response which leads to a time 
delay (z/c) for the discontinuities in the field, and the 
second part is the retarded response from the medium. 
In the following, we will prove that the second term is 
a continuous function. The discontinuity in the field is 
therefore determined by the first term, thus proving that 
the discontinuity always prapagates with the phase ve- 
locity. 

As the integral function in Eq. © is analytic in the 
upper half plane, J{£,;z) is equal to for ^ > since 



the integral function 
|w| — > oo. 
When ^ = 0, we have 
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where the integrations in the second and third lines are 
along the open semicircle C under the condition R ^ oo 
(see Fig. 1(a)). From the first to second lines, we have 
used that the intagration along the closed path composed 
of the real axis and the semicircle C (with R oo) is 
zero, because the integral, exp[—iz/c+ ^ "°"(^")^ ] _ jg 
analytical. In the third line we expanded exp[—iz/c + 
^ ujrj\uj)Z T^ — 1 for \uj\ OO, into ^ {Bn/uj^), where the 

^ n=l 

coefficients {-B„} are related to the coefficients {^n+i}. 



FIG. 1: (a) The integral path of Eq. ©; (b) The integral 
path of Eq. @. 



For ^ < 0, as shown in Fig. 1(b), we have 
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where the step leading to second line is because the func- 



tion 
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tends to zero when |a;| — s- oo; 

in the second line, the integral is along the closed path C 
(see Fig. 1(b)); in the third line the integrals are along 
the finite numbers of neighborhoods of the isolated sin- 
gular points and tangent lines, and these intergrals are 
finite; In the fourth line we describe a dashed circle of a 
radius r which embraces all these finite sigularities and 
tangent lines as shown in Fig. 1(b); and in the fifth line, 
we let ix) = r X e*'''. Therefore, we have the inequality 
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where M is the maximum value of the function 
^-i±^^iu,a{i^)z/c _ g«^?g«v on the dashed circle of 



3 



radius r. This means that J(^; z) is bounded. Because 
the function J{^,z) is bounded, the second term in Eq. 
© is continuous. A discontinuity in the input field am- 
plitude at time tg, will occur in the output field after 
passing through the medium only at the time equal to 
tg + z/c. This is the arrival time of information carried 
by the discontinuity after passing through the medium, 
which propagates with the phase velocity c. 

The intensity discontinuity might disappear under 
some specific cases due to the second term. The rea- 
son is that although |ap — \b\'^ ^ 0, |a-|-cp - |6-|-c|^ may 
be equal to zero for complex numbers a, 5, c. Once such 
a discontinuity disappears, there is no intensity disconti- 
nuity at any time, and the information is lost. It should 
be pointed out that there is no requirement for the form 
of the incident field in the above proof and particularity 
of the medium except the Kramers-Kronig relations. 

Next we consider the propagation of a nonanalytic dis- 
turbance in the derivatives of the field amplitude. Sup- 
pose the 0th to the {n — l)th order derivatives of the 
field are continuous functions, while the nth derivative 
has a discontinuity. How does this kind of nonalytical 
disturbance propagate? From Eq. JB)), we obtain 

&''E^iz,t) ,iAxd^Eo{z,t-z/c) , ,£^a"^,„(z,t) 
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we finally obtain 



a<"-^'5(i). Therefore, 
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where j;(e, z) = (-i)" dtje*"« {lu{- ■ ■ {uj{e^^ 



X e ^ — 1} — Bi} • • • } — Bn}- Following the same 
reasoning as discussed above, we can again prove that 
J„(^, z) is equal to for ^ > and is bounded for ^ < 0. 
Therefore, the second term in Eq. (|12() is a continu- 
ous function of time. The discontinuity in the nth order 
derivative of the output field solely comes from the first 
term in Eq. (fT^ (recall that the 0th to the (n — l)th 
order derivatives are continuous). 



FIG. 2: The intensities and their first order derivatives of 
the output fields after passing through the gain medium 
(solid lines) and the same distance of vacuum (dashed lines). 
The input field is E{0,t) = 2.0 exp[-tV(2cr?)] when t < 0, 
E(0,t) = 2.0 - 5.0t when < i < 0.2, and E{0,t) = 
exp[-(t - 0.2)^/(2(T?)] when t > 0.2. The parameter of the 
medium is the same as WKD's experiment jg. 



We can conclude that the discontinuity in the nth order 
derivative propagates with the phase velocity c. This dis- 
continuity might also be washed out under certain cases. 
Once this discontinuity disappears, there is no disconti- 
nuity at any other time (the information is lost). Any 
initial discontinuity in the nth order derivative (includ- 
ing the 0-th order) never leads to any new nonanalytical 
disturbances during the propagation. Therefore there is 
one to one correspondance between the disturbances (in- 
formation) in the output field and the input field. 

As a numerical example, the propagation time of the 
discontinuities in the first order derivative are calculated, 
see Fig. 2. There are two nonanalytical disturbances in 
the first order derivative of the input field at ti = 
and t2 = 0.2/is, and other parts of the input field are 
analytic. After passing through the gain medium of 
length L — 6cm, 8| these two discontinuities occur at 
time t[ 2 — ti,2 + L/c given by the phase velocity, but 
not by the group velocity. Please note that the initial in- 
tensity is a parabolic line in the interval (0, 0.2)/is. The 
information carried by the discontinuities propagate at 
the phase velocity c. 

In the second example, we calculate the propagation of 
the nonanalytical disturbances in the second derivative 
(see Fig. 3). There are two such discontinities at time 
ti = — 4.8^s and t2 — 4.8^s (see Fig. 3). The analyti- 
cal disturbances in the 2nd order derivative of the out- 
put field after passing through the gain medium arrive at 
t'l 2 ~ ti,2+L/c (same as passing through vacuum). That 
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— > oo). In this case, we can again show that the 
propagation of any nonanalytical disturbance is at the 
speed of c/uq (the phase velocity). 

We have proved that the information carried by the 
nonanalytical disturbances in the amplitude or in its any 
order derivatives (including the front of a pulse) propa- 
gates with the phase velocity, but not the group velocity. 
Our proof is based only on Maxwell's electromagnetic 
theory and the Kramers-Kronig relations without any 
condition on whether the medium is normal or anoma- 
lous. We believe this is why Einstein only considered the 
phase velocity in his relativity theory. 
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FIG. 3: The intensities and their 2nd order derivatives of 
the output field after passing through tire gain medium (solid 
lines) and the same distance vacuum (dashed lines). The 
input field is E{0,t) = A[l + cos(7rt/2.4)] when -7.2 < t < 
-4.8 and 4.8 < t < 7.2, E{0,t) = yl[1.5-|-0.5 cos(7rt/2.4)] when 
-4.8 <t< 4.8, and otherwise _E(0, t) = 0. The parameter of 
the medium is the same as WKD's experiment Ig. 

is to say, these two nonanalytical disturbances (the en- 
coded information) propagate through the gain medium 
with the phase velocity c. From the expression of the 
input field, we know that there are two more nonalytical 
disturbances in the third derivative of the input intensity 
at time ^3 = — 7.2/is and — 7.2/is. The disturbances in 
the third derivatives propagate also at the phase velocity 
c. 

In the two examples of numerical calculation, we used 
again medium of the type considered in the experiment 
[3 with anomalous dispersion. Here we would like to 
note that the Kramers-Kronig relations is applicable to 
the gain mediun. To our knowledge there is no dielectric 
media that do not obey the Kramers-Kronig relations. 
We also calculated the propagation time for the normal 
dispersive media. The calculations confirm the same re- 
sult that the nonanlytical disturbances propagate at the 
phase velocity, although the group velocities of these me- 
dia are smaller than the phase velocity. 

For a practical medium, the background refraction in- 
dex of the medium always exists. In this case, we can 
write n{uj) = no + n{uj) ( no is the backgroud refrac- 
tive index and is larger than unity, and n{uj) as 
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